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ABSTRACT 



We report the discovery of a new stellar companion in the KOI-13 system. KOI-13 is composed of two fast-rotating A-type stars of 
similar magnitude. One of these two stars hosts a transiting planet discovered by Kepler. We obtained new radial velocity measure- 
ments using the SOPHIE spectrograph at the Observatoire de Haute-Provence that reveal an additional companion in this system. 
This companion has a mass of between 0.4 M Q and 1 M Q and orbits one of the two main stars with a period of 65.831 ± 0.029 days 
and an eccentricity of 0.52 ± 0.02. The radial velocities of the two stars are derived using a model of two fast-rotating line profiles. 
From the residuals, we find a hint of the stellar variations seen in the Kepler light curve, which have an amplitude of about 1 .4 1 km s~ 1 
and a period close to the rotational period. This signal appears to be about three orders of magnitude larger than expected for stellar 
activity. From the analysis of the residuals, we also put a 3-cr upper limit on the mass of the transiting planet KOI-13.01 of 14.8 Mj up 
and 9.4 M Jup , depending on which star hosts the transit. We find that this new companion has no significant impact on the photometric 
determination of the mass of KOI-13.01 but is expected to affect precise infrared photometry. Finally, using dynamical simulations, 
we infer that the new companion is orbiting around KOI-13B and that the transiting planet candidate is expected to orbit KOI- 13 A. 
Thus, the transiting planet candidate KOI-13.01 is orbiting the main component of a hierarchical triple system. 

Key words, planetary systems - techniques: radial velocities - techniques: spectroscopic 



1. Introduction 

Among the 232 1 Kepler Objects of Interests (KOI) discovered so 
far (Batalha et al. 2012), KOI-13 is one of the most interesting. 
It is a common proper-motion binary with two bright (Va = 9.9 
mag and V# = 10.2 mag, Dommanget & Nys 1994) stars an- 
gularly separated by ~1.2" (Aitken 1904). One of these stars 
hosts a transiting planet candidate, KOI-13.01, orbiting every 
1.76 days (Batalha et al. 2012). A first analysis performed by 
Szabo et al. (2011) (hereafter Szll) determined that the plane- 
tary candidate is transiting the brightest component of the binary, 
KOI- 13 A. They also estimated the stellar parameters of both 
stars. They found that both the A and B components are A5-type 
stars with values of v sin/* of 65 ± 10km s _1 and 70+ lOkms -1 , 
respectively. From the analysis of the transit, they determined 
the radius of KOI-13.01 to be 2.2 ± 0.1Rj up and concluded that 
it is an inflated brown dwarf. They also found an asymmetry 
in the transit light curve that is characteristic of a misaligned 
companion transiting a hot and fast-rotating star affected by the 
gravity darkening effect (see Barnes et al. 2011, (Ball) and ref- 
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erences therein). By considering either beaming or ellipsoidal 
effect (Faigler & Mazeh 201 1), several studies have tried to mea- 
sure the mass of the transiting candidate using only the exquisite 
Kepler photometry. Shporer et al. (2011) (hereafter Shi 1) and 
Mazeh et al. (2012) (hereafter Mal2) found a mass of 9.2 + 1.1 
Mj up and 10 + 2 Mj up using the beaming effect, and Mislis & 
Hodgkin (2012) (hereafter M&H12)found a mass of 8.3 + 1.25 
Mj up using the ellipsoidal effect only. Thus, KOI-13.01 appears 
to be a massive and inflated hot jupiter orbiting one component 
of a binary system. In this letter, we present the results of our 
radial velocity follow-up observations of this system. 

2. SOPHIE observations 

2.1. Observations 

We performed spectroscopic observations of the KOI-13 system 
with the SOPHIE spectrograph (Perruchot et al. 2008; Bouchy 
et al. 2009) mounted on the 1.93-m telescope at Observatoire de 
Haute-Provence, France. We acquired 17 high resolution spectra 
from 201 1, March 25 to 2012, May 23 1 using the high efficiency 
mode (R ~ 39 000 at 550nm) of SOPHIE and the fast CCD read- 
out mode. The spectra were reduced with the online standard 
pipeline and the weighted cross-correlation functions (hereafter 
CCF) were computed using a F0 mask (Baranne et al. 1996; Pepe 
et al. 2002). 

1 prog. IDs: 1 1 A.PNP.MOUT, 1 1B.PNP.MOUT, 12A.PNP.MOUT 



1 



A. Santeme et al.: An additional companion in the KOI- 13 system 



2.2. Cross-correlation function modeling 

Since the two main stars of the KOI- 13 system are separated by 
~ 1.2", their fluxes were observed simultaneously within the 3" 
diameter of the SOPHIE fiber. Thus, the spectra of both stars are 
blended, and there is no straightforward method for identifying 
the two contributions. To measure the radial velocity (hereafter 
RV) of both stars in the KOT13 system, we first normalized the 
CCF from the blaze response of the spectrograph. The normal- 
ized CCFs are plotted in Fig. 1. We then computed a rotational 
profile to fit the data. The observed profile of a rapidly rotat- 
ing star is the result of the convolution of a rotational profile 
G{v) with the instrumental profile (Gray 2005). For the high ef- 
ficiency mode of SOPHIE, we used a Gaussian with a FWHM of 
10km s as the instrumental profile. The rotational profile G(v) 
was defined to be (Gray 2005, p. 465) 



G(u) = 



2 (1 - e) Vl - (v/ul) 2 + \ns [l - (v/v L ) 2 ] 



nv L {\ - e/3) 



(1) 



where Vl is identical to the sky-projected rotational velocity 
v sin i+ and s is a limb-darkening coefficient. For the CCF mod- 
eling, we fixed the limb-darkening coefficient to e = 0.3 (see 
Gray 2005, Fig. 17.6). 

2.3. Distinguishing the radial velocities of both components 

To adjust the CCFs, we thus constructed a model with two fast- 
rotating stars. We fitted the 17 CCFs simultaneously with this 
model of two fast-rotating stars allowing both the RV, contrasts 
and L>sin/* to vary. We adjusted the contrast of the two stars 
for each CCF independently since both the seeing and the tele- 
scope guiding may favor one or the other star during the expo- 
sure. The best-fit model of each CCF as well as the residuals 
are displayed in Fig. 1. One of the two stars seems to have a 
CCF profile with a larger area. We know from Hilditch & King 
(1986) that the relative profile area of a binary can be used to 
measure its luminosity ratio. However, this is not enough evi- 
dence to be able to match each profile with the known star. In 
consequence, we refer to the component with a slightly larger 
(smaller) profile area as KOI- Off (KOI-13yS, respectively) to 
avoid any confusion. We find i/sin/* a = 76.6 + 0.2kms _1 and 
i/sin/* a = 62.7 + 0.2 km s -1 . We note that these values agree 
with the values found by Szl 1 to within l-cr. The median values 
for the contrasts are c a = 0.58 + 0.09% and c = 0.59 ± 0.05%. 

The derived RV and their uncertainties for both stars are 
listed in Table 1 and displayed in Fig. 2. We found no clear RV 
variation for the star KOI-13ff with a systemic RV of 1.20 + 
1.24kms~ 1 . On the other hand, KOI-13/3 has a large amplitude 
RV variation. We performed a Keplerian fit to the data of KOI- 
13/3 using a Levenberg-Marquardt algorithm. The best solution 
has an eccentricity of 0.52 ± 0.02, a period of 65.831 ± 0.029 
days, and a semi-amplitude of 12.42 + 0.42 kms~'. The param- 
eters are listed in Table 2 with their uncertainties. This variation 
is caused by the reflex motion of the star orbited by a stellar 
companion with a minimum mass of 0.40 ± 0.03 M . We called 
this new stellar companion in the system KOI-13y. We estimated 
a detection limit for a third component in the CCFs, after re- 
moving the main contributions of both A5 stars. The full set of 
SOPHIE spectra could exclude the presence of a third Gaussian 
non-rotating profile with an amplitude of 0.4%, with more than 
10-cr confidence. Within the first-order assumption that the flux 
ratio scales as the CCF area, one can derive an upper-limit to 
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Fig. 1. Upper panel: Individual CCFs superimposed on a two- 
star CCF model (black lines). The CCFs are arbitrary shifted 
down with increasing time. Lower panel: The CCF residuals dis- 
played in the same way as in the upper panel. The two vertical 
lines indicate the limits of the CCFs. 



Table 2. KOI-13y orbital parameters 



Orbital period P [d] 


65.831 i 


: 0.029 


Epoch of periastronT p [BJD-245000] 


5943.155 


±0.181 


Eccentricity e 


0.52 + 


0.02 


Argument of periastron co [°] 


32 i 


:3 


Systemic radial velocity u [kms -1 ] 


3.98 ± 


0.24 


Radial velocity semi-amplitude K [kms -1 ] 


12.42 d 


:0.42 



the mass of KOI-13y of 1 M G . With less conservative assump- 
tions, this upper limit decreases to 0.7 M G - when 1 1 of the 17 
residual CCFs can exclude a third component at the 3-cr limit. 
In the following, we choose the conservative approach and an 
upper limit of 1 M Q . KOI-13y is not transiting its host star since 
a deep transit would have been seen in the Kepler light curve. 
Considering the whole possible mass range for KOI-13y and the 
expected radii (Carroll & Ostlie 1996), we estimated that an or- 
bital inclination of up to 87° is compatible with the absence of 
a transit (eq. 7, Winn 2010). Using eq. 1 from Faigler & Mazeh 
(2011), we estimated that KOI-13y should produce a beaming 
effect on KOI-13/? at the level of about 165ppm, without any di- 
lution. This amplitude may be detectable with a careful filtering 
of the Kepler light curve that is beyond the scope of this letter. 
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Fig. 2. Radial velocity curves and residuals of KOI- 13a (blue squares) and KOI-/? (red circles) measured by SOPHIE against time 
(left panels) and phase-folded to the best-fit period of KOI-13/? (right panels). The dashed blue line represents the systemic RV of 
KOI- 13a. 

2.4. Stellar variability 

Another signal in the Kepler light curve varies with a semi- 
amplitude of about ~20 ppm and a period of 1.059 + 0.002 d 
(Mal2; Shi 1 ; Szabo et al. 2012). The discoveres of this signal 
claimed that the signal has a stellar origin. Szabo et al. (2012) 
also claimed that the signal has a period close to the expected 
rotational period of the star and may be either granulation noise 
or starspot features (see also Balona 201 1). 

We phase-folded the residuals of both stars to the expected 
stellar-rotational period within the errors. We found that the 
residuals of KOI- 13a appear to have a circular variation with a 
semi-amplitude of 1.41 +0.38 kms -1 (3.6-<t significance level). 
The best-fit period of this signal detected by RV is 1.0642 + 
0.0003 d, at 2.6-cr from the photometric period. This amplitude 
is slightly larger than the expected variation caused by the planet 
candidate of about 1 tans" 1 . The residuals of KOI- Do-, phase- 
folded at the stellar rotational period are plotted in Fig. 3. No 
similar signal is seen in the residuals of KOI-13/?. If the varia- 
tion seen in the light curve is due to spots that produce a depth 
of about 60ppm, we might have expected to measure a RV am- 
plitude at the level of a few ms _1 (Aigrain et al. 2012). If the 
signal seen in RV is real, it should have a different origin from 
starspots and might represent stellar pulsation (Shi 1), as in the 
case of WASP-33 (Collier Cameron et al. 2010). 




Fig. 3. Radial velocity residuals of KOI- 13a phase-folded to the 
best-solution period (P = 1.0642 d). The black line is the best-fit 
circular model that fits the data. 



Fig. 4. Radial velocity residuals of KOT13a (left panel) and 
KOI-13/? (right panel) phase-folded to the ephemeris of the 
Kepler transiting-planet candidate KOI-13.01 (P~ 1.76 d). The 
solid, dashed, and dot-dashed lines represent the 1-cr, 2-cr, and 
3-cr (respectively) contours of a circular orbit. 

3. On the mass of the Kepler transiting candidate 
KOI-13.01 

3.1. Radial velocity upper-limit 

We used the residuals of the two components a and /? to place 
constrains on the mass of KOI-13.01 detected by Kepler. We 
phase-folded the residuals (see Fig. 4) at the Kepler ephemeris 
(P=l. 763589 d, t =54953.56498 BJD). The 3-cr upper limit is 
14.8 Mj up or 9.4 M Jup , depending on whether the host star is 
KOI- 1 3c or KOI-13/?, respectively. We note that to be conser- 
vative, we did not remove the stellar variability from KOI- 13a. 
Removing this signal, we found a slightly lower upper-limit to 
the mass of KOI-13.01 of 14.4 Mj up . The asymmetric transits of 
KOI-13.01 shown by Sz 11 and Ball must be produced by a hot 
and fast rotating star, completely ruling out KOI-13y as the host 
of the transiting planet candidate. If KOI-13.01 were transiting 
KOI- 13a, our results would agree with the mass derived using 
only photometry by Shi 1, Mal2, and M&H12 and would con- 
firm the planetary nature of KOI-13.01. On the basis of an upper 
limit of 9.4Mj up at 3-cr compared to the photometric determi- 
nations, the scenario in which KOI-13.01 is transiting KOI-13/? 
seems more unlikely. 



3.2. Impact of KOl-13y on the photometric mass 
determination of KOI-13.01 

The mass estimate of KOI-13.01 was performed using the am- 
plitude of the beaming or ellipsoidal effects (Shi 1; Mal2; 
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M&H12). This photometric determination does not take into ac- 
count the amount of light from KOI- 1 3y. Equations 6 and 7 of 
the beaming and the ellipsoidal effects provided by Mal2 show 
that the mass is proportional to the amplitude of these effects. 
The flux from KOTOy impacts linearly on the mass determina- 
tion of KOI- 13.01. 

We investigated the impact of this new stellar companion 
on the mass of the planetary companion. We computed the di- 
lution factor that affects the light curve of KOI- 13 A, which is 
expected to be the host of the transiting planet (Szll). We gen- 
erated Planck curves of two A-type stars with radii of 2.55 R 
and 2.38 R G and T eff of 85 1 1 K and 8222 K, respectively (Szll). 
We generated a third Planck curve for KOTOy, assuming a true 
mass, radius, and T e ff according to table B.l of Carroll & Ostlie 
(1996). We integrated the three Planck curves within the Kepler 
and Spitzer IRAC at 3.6jum and 4.5/im bandpasses. We allowed 
the mass of KOTOy to vary between 0.4 M and 1 M Q as its 
inclination is unknown, but constrained by SOPHIE spectra. We 
found that KOTOy increases the dilution factor affecting KOT 
13A by less than 0.8% within the Kepler bandpass. Thus, the 
photometric measurement of KOT13.01's mass is insignificantly 
affected. However, in the infrared, the additional dilution factor 
produced by KOTOy reaches up to 2.2% for both sets of IRAC 
dara at 3.6/mi and 4.5pm. Thus, this third star is expected to af- 
fect the precise characterization of the planetary atmosphere. 



4. Constraining the scenarios 

With an angular separation of 1.2", the two main components 
of the KOTO system were observed simultaneously within the 
3" diameter fiber of the SOPHIE spectrograph. We detected a 
new companion in the system, called KOTOy. We were un- 
able to measure the mass of the transiting-planet candidate KOT 
13.01 nor determine about which star, a or B, it orbits. Szabo 
et al. (2011) determined that the transiting planet candidate 
KOTO. 01 is transiting the star A using a ground-based lucky- 
imaging photometry. This suggests two different scenarios: 1) 
both companions orbit the same star, i.e. B = A and a = B 
or, 2) both companions orbit a different star, i.e. a = A and 
B = B. From the line profile area, we assumed that KOTOa- 
is the brightest component. The v sin i+ values derived from the 
CCF fit also hints that a = A and B = B. 

To more tightly constrain the two scenarios, we performed 
dynamical simulations using the Mercury6 code (Chambers 
1999), assuming that KOI-13.01 and KOI-Oy orbit the same 
star (scenario 1) and allowing different inclinations and masses 
for KOI- 1 3y. As the inclination of the orbital plane of KOI- 1 3 .0 1 
(85.9°, Ball) is of the same order as the maximum inclination 
of KOTOy (87°, see section 2.3), the orbits can be coplanar. 
In this configuration, the eccentricity of KOTO. 01 can be ex- 
cited between and 0. 15 with periodical variations over about 30 
years. The time of the occultation found in the Kepler light curve 
(Szll) indicates that the orbit is highly circular (Mal2). From 
our simulations, we found that a non-zero eccentricity will be 
detected within the lifetime of Kepler. However, we also found 
that such a configuration would produce transit-time variations 
(TTV) at the level of several minutes. We analyzed the Q2 and 
Q3 short-cadence Kepler light curve and found no TTV at the 
level of 24 seconds, at 3-cr, using a procedure similar to Barros 
et al. (2011). If KOI-Oy has a different inclination than KOI- 
13.01, the orbit of the latter is either unstable (for a mass of KOI- 
Oy > 0.6 M ) or should display significant TTV, non-zero ec- 
centricity, as well as inclination variations on a short timescale. 



We thus conclude that KOI-Oy and KOTO. 01 are orbiting two 
different stars, which confirms that a = A and B = B. 

5. Conclusion and discussions 

We have performed RV follow-up observations of the KOTO 
system with SOPHIE. The KOTO system is dominated by two 
bright A-type stars with v sin/* of approximately 70km s -1 . By 
modeling the observed CCFs using two rotational profiles, we 
have been able to distinguish the RVs of the two stars and found 
that this system hosts an additional companion in the stellar 
regime, KOI-Oy. This companion has a mass of between 0.4 
M Q and 1 M , an orbital period of 65.831 + 0.029 days and 
an eccentricity of 0.52 ± 0.02. This previously unknown stel- 
lar companion does not significantly increase the dilution of the 
Kepler light curve, thus does not affect the photometric mass de- 
termination of KOI- 13.01. However, this companion is expected 
to make a significant contribution to infrared bandpasses. In the 
residuals of KOTOff, we found a hint of stellar variability with a 
period of about 1.064d, which is close to the period determined 
by photometry. This variation has an amplitude of 1.41 ± 0.38 
kms and is slightly longer than the expected variation due to 
the planet candidate KOI-13.01 (K < 1 kms 1 ). If this signal is 
real, constraining the mass of KOI-13.01 using RV data will be 
challenging. We have also provided a 3-cr upper-limit constraint 
based on the RV residuals of about 14.8 Mj up or 9.4 Mj up on the 
mass of KOI-13.01, depending on whether KOI- Off or KOI- 13/? 
is transited, respectively. These limits in mass are compatible 
with the photometric measurements performed by Shi 1, Mal2, 
and M&H12. KOI-13.01 thus appears to be a massive transiting 
planet orbiting a member of a triple system. 

Finally, we used dynamical simulations of the KOTO sys- 
tem with the Mercury 6 code (Chambers 1999) to show that both 
companions cannot orbit the same star and therefore conclude 
that KOI-Off is KOI-OA and KOI-OyS is KOI-OB. 
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